SHARP ASYMPTOTICS OF METASTABLE TRANSITION TIMES FOR 

ONE DIMENSIONAL SPDES 



FLORENT BARRET 



Abstract. We consider a class of parabolic semi-linear stochastic partial differential 
equations driven by space-time white noise on a compact space interval. Our aim is to 
obtain precise asymptotics of the transition times between metastable states. A version 
of the so-called Eyring-Kramers Formula is proven in an infinite dimensional setting. 
The proof is based on a spatial finite difference discretization of the stochastic partial 
differential equation. The expected transition time is computed for the finite dimensional 
approximation and controlled uniformly in the dimension. 



1. Introduction 

Metastability is a phenomenon which concerns systems with several stable states. Due 
to perturbations (either deterministic or stochastic) the system undergoes a shift of regime 
and reaches a new stable state (see e.g. [15] by Cassandro, Galves, Olivieri and Vares, the 
book |36j by Olivieri and Vares and the lecture notes [6] by Bovier). Typical examples of 
metastable behavior can be found in chemistry, physics (for models of phase transition) 
and ecology. 

In this article, our aim is to understand metastability for a class of stochastic partial 
differential equations. We consider the Allen-Cahn (or Ginzburg-Landau) model which 
represents the behavior of an elastic string in a viscous stochastic environment submitted 
to a potential (see e.g. Funaki [27]). This model has other interpretations in quantum field 
theory (see [211 HB] and the references therein) and in statistical mechanics as a reaction 
diffusion equation modeling phase transitions and evolution of interfaces (see Brassesco 
and Butta [I21[l3]). 

More precisely, we deal with the following equation, for (x, t) G [0, 1] x M"*" 



dtu{x, t) = 79^^m(x, t) - V'{u{x, t)) + V2eW (1.1) 

where 7 > 0. TV is a space-time white noise on [0, 1] x in the sense of Walsh [IQ] and 
£ > is the intensity of the noise. V is a. smooth real valued function on M called a local 
potential. We consider two boundary conditions: Dirichlet boundary conditions (for all 
t e M"*", u{0,t) = u{l,t) = 0) and Neumann boundary conditions ((9^.u(0,t) = dxu{l,t) = 
0). The initial condition is given by a continuous function uq which satisfies the given 
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boundary conditions. Existence and uniqueness of an Holder-continuous solution in the 
mild sense have been proved by Gyongy and Pardoux in [29] . 

Paris and Jona-Lasinio in [2T] are among the first ones to analyze Equation (II. ip for a 
double well potential 

4 2 

In this case, V has only two minima which are +1 and —1. One expects that the model (II. ip 
has several stable states and that a metastable behavior occurs. The authors introduced 
a functional potential S and interpreted (11.11) as the stochastic perturbation of an infinite 
dimensional gradient system: 

dtu = -— + V2^W (1.3) 
where for a differentiable function, 

S{<P) = l^'^\<P'{x)\' + V{<P{x))dx. (1.4) 

5* represents the free energy. |^ is the Prechet derivative of S i.e. the infinite dimensional 
gradient of S. 

For more general functions V (real valued functions), we can define a similar potential 
S as in (ll.4p which determines a potential landscape. Under the stochastic perturbation, 
this potential landscape is explored by the process u defined in (11.11) . While the system 
without noise (i.e. e = 0) has several stable fixpoints (which are the minima of S), for 
e > transitions between these fixpoints will occur at a suitable timescale. The transition 
paths go through the lowest saddle points. Thus, minima and saddle points of 5* have a 
key role to understand metastability but it is often a hard task, given a potential V (and 
thus S), to completely compute and comprehend the geometrical structure of the energy 
landscape. However, some elegant method exists (see e.g. [221 HTj ). 

The model (ll.3p is an infinite dimensional generalization of the finite dimensional systems 
investigated by Preidlin and Wentzell [25] and by Bovier, Eckhoff, Gayrard and Klein 
in [21 [To]. Moreover, we will see that (ll.ip is rigorously the limit of a gradient finite 
dimensional system (via a spatial finite difference approximation). 

Our aim is to derive precise asymptotics of the expected transition time i.e. the time 
needed, starting from a minimum (pQ of S, to hit a set of lower minima. We define the 
hitting time Ts{B) by t^^B) = inf {t > 0, u{t) G B} where S is a disjoint union of small ball 
around some minima of S lower than 0o- We prove that the expected time, K^^[ts{B)], has 
a very distinctive form known as the Arrhenius equation (Theorem 12. 6p . This expectation 

E^,[T,{B)]=Ae^/^{l + 0{V^\He)f')) (e ^ 0) (1.5) 

where E is the activation energy and A is the prefactor. E has been computed by Paris and 
Jona-Lasinio for the double well potential (11.21) using a large deviation approach (Theorem 
1.1 [21]). E is exactly the minimum height of potential that a pathway has to overcome 
to reach B starting from (pQ. The prefactor A is a constant (for our set of hypotheses) 
and depends only on the local geometry of the potential S near the minimum 0o and near 
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the passes (or saddle points) from 0o to the set B. The order ©(a/^ |ln(£)| ) of the error 
term comes directly from the local approximation of the potential 5* by its quadratic part. 

For the double well potential (11. 2p with Neumann boundary conditions, Paris and Jona- 
Lasinio proved that S has only two global minima, denoted m and — m (correspond- 
ing roughly to the constant functions 1 and —1 resp.). For some 7, this model has a 
unique saddle point a = (the constant function 0). We deduce from Theorem 12.61 that 



K_m[Ts{B~^)], for a small ball B"^ in the suitable norm around m, takes the form (11. 5p with 
E = S{a) — S{—m) and 



where {Xk{((>))k^i are the eigenvalues of the second Frechet derivative of the potential S 
at a point and A~((t) is the unique negative eigenvalue at the saddle point a. Using 
asymptotic expansion of the eigenvalues, we prove that the infinite product converges. It 
is exactly the equivalent for an operator of the classical determinant of a matrix. We also 
mention the fact that this infinite product has a nice expression in terms of solutions of 
linear differential equations (see e.g. Levit and Smilansky [32]). 

Eyring in [20] and particularly Kramers in investigate the case of a one dimensional 
diffusion as a model for chemical reaction rates and express rates instead of expectations. 
Their formula is known as the Eyring- Kramers Formula. It takes the form (II. 5p with the 
prefactor given by a formula similar to (II. 6p but with a single factor in the product (there 
is only one eigenvalue). 

Similar Eyring- Kramers Formulas exist through a wide range of reversible Markovian 
models from Markov chains, stochastic differential equations. For finite dimensional dif- 
fusions, Freidlin and Wentzell in [23], proving that these systems obey a large deviation 
principle, obtained the activation energy in terms of the rate function. In recent years, the 
potential theory approach initiated by Bovier, Eckhoff, Gayrard and Klein in [HI [ID] has 
allowed to give very precise results and led to a proof of the Eyring-Kramers Formula for 
gradient drift diffusions in finite dimension. Moreover, the potential approach originate 
from Markov chains (see [ZJ El [6]) and have been refined to obtain metastable transition 
times for specific models (see e.g. [5| fTTj). 

Formula (II. 6p is then the extension of the Eyring-Kramers Formula to a class of one- 
dimensional SPDEs (ll.ip . Maier and Stein in [33] obtained heuristically this formula and 
Vanden-Eijnden and Westdickenberg in [32] used it to compute nucleation probability. 

Specifically, the system (II. ip and its metastable behavior have been studied for at least 
thirty years using mainly large deviation principle and comparison estimates between the 
deterministic process ( (ll.ip with e = 0) and the stochastic process defined by (II. ip . Cas- 
sandro, Olivieri, Picco [16] obtained similar asymptotics as Faris and Jona-Lasinio [21] 
when the size of the space interval is not fixed and goes to infinity as e goes to suffi- 
ciently slowly. These results first prove the existence of a suitable exponential timescale 
in which the process undergoes a transition. 



A 



n 



Afe(a)| 



(1.6) 
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In the same case as fll.2p . Martinelli, Olivieri and Scoppola [31] obtain the asymptotic 
exponentiahty of the transition times (Theorem 4.1 [21]). Also, Brassesco [12] proves that 
the trajectories of this system exhibit characteristics of a metastable behavior: the escape 
from the basin of attraction of the minimum — m occurs through the lowest saddle points 
(Theorem 2.1 [T2|) and the process starting from — m spends most of its time before the 
transition near — m (Theorem 2.2 |12]). 

In this paper, we consider a local potential V (satisfying Assumptions 12. ll and 12. 4( ) and 
we rigorously prove an infinite dimensional version of the Eyring-Kramers Formula. Our 
method relies on a spatial finite difference approximation of Equation (11. ip introduced by 
Berglund, Fernandez and Gentz in [3l H] as a model of coupled particles submitted to a 
potential. The computation of the expected transition time for the approximated system 
gives us the prefactor, the activation energy and some error terms. We need to control 
the behavior of these error terms as the step of discretization goes to (or equivalently as 
the dimension N of the approximated system goes to +oo). To this aim, we adapt results 
from |2] by Bovier, Meleard and the author. 

As proved by Funaki [27] and Gyongy [28], the solution of the approximated system 
converges to the solution of the SPDE. By combining different results from SPDE theory, 
large deviation theory (from Chenal and Millet [17]) and Sturm-Liouville theory we are 
able to take the limit of the finite dimensional model in order to retrieve the SPDE (11. ip . 
We also need to adapt estimates on the loss of the memory of the initial condition (from 
Martinelli, Scoppola and Sbano [SUES]) uniformly in the dimension. 

The use of spatial finite difference approximation is quite natural since we consider our 
SPDEs in the sense of Walsh [10], limited to the case of space-time white noise. Other 
approximations could be possible, notably the Galerkin approximation should lead to 
similar results for a different class of SPDEs in the framework of Da Prato and Zabczyk 
(see the book [T9]). 

The article is organized as follows. In Section |2l we present the equation, the assump- 
tions, the main theorem (Theorem 12. 6p and a sketch of its proof. Then in Section |3l we 
adapt the convergence of the approximations and prove convergence of the approximated 
transition times. In Section m we state large deviations estimates by Chenal and Millet 
[T7] . contraction results by Martinelli, Olivieri, Scoppola and Sbano [Ml |25] and prove 
a uniform control in the initial condition uniformly in the dimension. In Section El we 
recall results about eigenvalues and eigenvectors of Sturm-Liouville problems and prove 
the convergence of the prefactor. In the last section, we compute the expected transition 
times uniformly in the dimension. 

We will use the following notations henceforth. For a functional space C, equipped with 
a norm \\-\\(^, we denote by Che the closed subspace in the C topology of the functions in C 
satisfying the suitable boundary conditions (Dirichlet or Neumann). For / G L°°([0, 1] x 
[0,T]) we set the norm of this space ll/llo^y or simply \\f\\^ when T = +oo. 

Acknowledgments. I am very grateful to Anton Bovier and Sylvie Meleard for suggesting 
this topic and for constant help and advice. I am indebted to the Hausdorff Center for 
Mathematics Bonn for financial support of numerous visits to Bonn. Part of the work 
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2. Results 

2.1. The Equation. The assumptions are of two kinds: some on the local potential V, 
others on the functional potential S. We first start with the hypotheses on V. 

Assumptions 2.1. We suppose that: 

• V isC^ on R. 

• V is convex at infinity: there exist R,c> such that for \u\ > R 

V"{u)>c>0. (2.1) 

• V grows at infinity at most polynomially: there exist p,C > such that 

V{u) < C{1 + \uf). (2.2) 

These hypotheses are made to avoid complications for the definition of the solution u of 
(11. ip and to allow the computations of the derivatives of S. 

Let {Q, P) be a probability space on which we define a space-time white noise W as 
defined in |1D] equipped with a filtration {J^t)t^o with the usual properties. The integrable 
processes for the white noise are the predictable measurable processes in L^(r2xM_|. x [0, 1]). 
We denote by gt{x,y) the density of the semi-group generated by ■ydxx on [0, 1] with the 
suitable boundary conditions. 

Let us recall that a random field m is a mild solution of (II. ip if 

(1) -u is almost surely continuous on [0, 1] x and predictable 

(2) for all {x,t) G [0, 1] x M+ 

u{x,t) = / gt{x,y)uo{y)dy - / / gt^,{x,y)V'{u{y, s))dyds 



JO JO JO 

+ V2^ [ [ gt.s{x,y)W{dy,ds). (2.3) 
Jo Jo 

We state from [29] the following result on the existence, uniqueness and regularity of 
the solution. 

Proposition 2.2 ([29]). For every initial condition uq G Cfec([0, 1]), the stochastic partial 
differential equation (II. ip has a unique mild solution. Moreover for all T > and p ^ 1, 



E 



sup \u{x,t)\^ 

[0,T]x[0,l] 



^C{T,p). (2.4) 



The random field u is essentially ^-Holder in space and ^-Holder in time. 

The only complication comes from the fact that V is not globally Lipschitz but prevents 
the process to go to infinity. From Assumptions 12.11 we have 

- xV'{x) < C. (2.5) 
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The proof of Proposition 12.21 is standard and uses mainly estimates on the density gt{x, y). 

Remark 1 . The definition of the stochastic convolution (the last expression of the right- 
hand side of (12. 3p ) requires the density of the semi-group to be in L^([0, 1] x [0,T]) for 
every T > 0. Unfortunately, that is only true in dimension one. For higher dimensions, the 
stochastic convolution does not define a classical function but a distribution in a Sobolev 
space of negative index [ID] . 



2.2. Stationary Points. As for the finite dimensional case, the minima and saddle points 
of S play a crucial role. To this end, we first specify what is the "gradient" (or the Frechet 
derivative) of the functional S. Let us recall that S is defined, for G H^^, by 

S{<P) = l^'^\<P'{x)\' + V{<P{x))dx. (2.6) 
For (j), h in C|^([0, 1]) we have a Taylor expansion of S at the second order in h 

S{<f> + h) = 5(0) + D^S{h) + ^DlSih, h) + 0{\\h\\%) (2.7) 

where 11/^11^-2 = ll^lloo"'"ll^'lloo"'"ll^"lloo- By integration by parts we compute the differentials 
D^S and D'^S. The first order differential is a linear functional which takes the form 

D^S{h) = [ [-70" (x) + V'{(P{x))]h{x)dx. (2.8) 

The Frechet derivative is = — 70"(x) + V'{(f){x)). The second order derivative (the 
Hessian operator) takes the form 

DlS{h,h)= [ h{x)[--fh"{x) + V"{(P{x))h{x)]dx. (2.9) 
Jo 

We denote by H^S* the Hessian operator at 0: 

n^Sh{x) = -jh"{x) + V"{(P{x))h{x). (2.10) 

The Hessian operator is a Sturm-Liouville operator. 

We say that is a stationary point of S" if is solution of the non-linear differential 
equation 

^ = -70" + V^'(0) = 0. (2.11) 
Let us now fix two points 0, G Cbc{[0, 1]) and define some quantities. 

r(0 ^ v^) = {/, /(o) = 0, /(I) = ijje C([o, 1], CU[0, i]))} (2.12) 

is the set of continuous paths from to tp. For / G r(0 tp), f denotes the set of maxima 
of the path /, 

/ = I fito),to G argmax5(/(t)) I . (2.13) 

[ t6[0,l] J 

The saddle points are passes from a valley to another one. The definition uses this idea. 
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Definition 2.3 (Saddles). For any G Cbc([0,l]), we define 5(0, the minimum 
heiglit needed to go from (p to 

5(0, ij) = S{tP, 0) = inf {5(0), G /, / G r(0 ^ V')} • (2.14) 

For 0,-0 sucli that S{(j),i/j) < oo, we denote ^(0,-0) the set of admissible saddles: the 
points which realize the maximum along a minimal pathway 

^(0, i^) = {ae CU[0, 1]), S{a) = 5(0, z^), 3/ G r(0 -> ^), a G /} . (2.15) 

The set of admissible saddle points is very important to compute the prefactor of the 
mean transition times. Near these points the process spends the most crucial time as it 
passes from a basin of attraction to another one. 

We now present the assumptions on 5. 

Assumptions 2.4. We suppose that: 

• 5 has a finite number of minima and saddle points. 

• All the minima and saddle points of 5 are non-degenerate (i.e. hyperbolic): at each 
point, the Hessian operator has non-zero eigenvalues. 

Assumptions 12.41 are structural. The finite number of stationary points provides a sim- 
ple generalization of the case where there is only one saddle point. The non- degeneracy 
condition is necessary in order to approximate locally at the minima and saddle points 
the potential by its quadratic part. If this is not the case the prefactor in fll.Sp is not a 
constant but should have a dependence in e. 

Connections between Assumptions 12.11 and 12.41 are not straightforward. Proving that a 
given potential 5 satisfies Assumption 12.41 is not easy, a precise analysis is often needed. 
Moreover if we want to investigate the dependence of the potential 5 on the parameter 7, 
bifurcations can occur and the landscape do not satisfy Assumption 12.41 for some critical 
values of 7. See Berglund, Fernandez and Gentz [21 H] for the finite and infinite dimensional 
cases for the double well potential. However, results exist (see [1] and references therein) 
on the generality of Assumption 12.41 

In addition, under Assumptions 12.41 and 12. 1^ the deterministic dynamical system (i.e. 
f ll.ip without the white noise) satisfies a Morse-Smale structure (see [T^ |22] and the 
references therein). This means that the attractor of the dynamical system consists of 
equilibria and heteroclinic orbits connecting these equilibria. Methods has been developed 
by Fiedler and Rocha in [22], by Wolfrum in [41] to compute the global attractor of the 
deterministic system. 

Remark 2. is the convenient functional space for the process since 5(0) < +00 if and 
only if is in H^{[0, 1]). In fact from the upper bound (12. 2 p and lower bound (12. ip on V 
we get 

Cimim - 1) ^ Si<P) ^ C'MWm + UWm + !)• (2-16) 
Each function in i?^([0, 1]) is continuous and even a-Holder continuous (for < a < |). 

For each G C([0, 1]), we define the quantity Det{'H^S): 
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• for Dirichlet boundary conditions, let / be the solution on [0, 1] of 

'H^Sf = /(0) = 1 /'(0) = (2.17) 

then Det{H^S) = /(I) 

• for Neumann boundary conditions, let / be the solution on [0, 1] of 

'H^Sf = /(0)=0 /'(0) = 1 (2.18) 

then BetiH^S) = /'(I). 
Let us recall that, as a regular Sturm-Liouville operator, H^S has a countable number 
of eigenvalues, all of them real. We denote by {Xk{(p))k^i the sequence of these eigenvalues 
in the increasing order. The definition of DetlH^S) is justified by the following lemma. 

Lemma 2.5 ([32]). For any and ip with non- degenerate Hessian operator, the infinite 
product Y[T=i Afc(^) convergent and we have 

n Afc(0) _ Det(-H^^) 



This lemma relates the infinite product of the ratio of eigenvalues to a ratio of terminal 
values of solutions. We find an elementary proof in [32] by Levit and Smilansky which relies 
on two different expressions of the Green function associated to the problem T-L^Sf = 
satisfying the boundary conditions. In fact, the Green function could either be expressed 
using the spectral decomposition of HfpS or expressed as a linear combination of two well- 
chosen fundamental solutions (of the second order linear differential equation). 

2.3. Main results. Before stating the main result, we describe the set of minima and 
saddle points. In fact, the prefactor depends greatly on the geometry of a graph connecting 
the minima to each other through the saddle points (so-called the 1-skeleton connection 
graph by Fiedler and Rocha in [23]). We define this graph and express the prefactor partly 
as an equivalent conductance on this graph. 

We denote by M the set of minima of S. Since by Assumption 12. 4^ there is a finite 
number of stationary points, we order the minima by increasing energy. We denote by 
(pi, (1)2, . . . , (pm, 'nT' = \-M.\, the different minima indexed by increasing energy 

^(0l)^^(02)^...^^(0j. (2.20) 

We denote by Aii, the subset of minima A^/ = 02, • • • , <t>i} for 1 ^ / ^ m. 

We consider the transitions from a minimum to Aii for I < Iq. These are the only 
visible metastable transitions. We will see from large deviations estimates, that to go from 

a minimum </> to another ip, it requires a time of order exp {^S{(t),il)) — S{(j))/e^. The time 

required to make the reverse transition is also of order exp {^{i), (f>) — S{ip)/e^ . Therefore 
if S{ij) > S{(j)), we get 

^(0,^)-^(0)>^(^,0)-^(^) (2.21) 

and the time required to go from to ^/^ is much larger than for the reverse transition. So 
we cannot see the reverse transitions since there are absorbed by the direct ones. If some 
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minima have the same potential, we can suitably order them to consider a transition from 
one minimum to another one at a same height. 

Let us now construct the weighted graph of paths from to A^^. We denote S = 
S{(j)iQ,Aii) the common potential of the saddles. The minima Ai are the vertices of the 
graph, the saddle points in ^(^/p, A^/) are the edges. We connect an edge a between two 
vertices G if the saddle a is a pass between the valleys of (p and ip: there exists 
/ G r(0 — )■ -i/;) such that / has a unique element and f = a. Existence of this graph is 
ensured by Assumptions 12.41 (see [23] and references therein). 

Each saddle point in ^(0;^, A^;) has a unique negative eigenvalue from the Morse-Smale 
property and the hyperbolicity of the stationary points. The weight associated to an edge 
a is defined as 

wia) = (2.22) 

where A~((t) is the unique negative eigenvalue of H^S. 

and a~ denote the two minima connected by a given edge a. Let us recall that we 
have m minima in Al. For a real valued vector a G indexed by the minima in Al, we 
consider the following quadratic form 

Q{a)= J2 w{a){a{a+)-a{a-))^. (2.23) 

We define C*{(f>ig,J^i) the equivalent conductance of the graph between 0;^ and Al; as 

C*{(f)i,,Mi) = inf {Q(a), a G M'", a(0,J = 1, a(</.) = 0, G Al^} . (2.24) 

This conductance is an approximation of the capacity between a neighborhood of (pi^ and 
Al/. In some sense, we replace the continuous landscape defined by 5 by a graph containing 
the relevant geometric structure of the landscape. 

Let us denote by Bp{(f)), for </> G H^cl'^, 1], the ball of center and radius p in H^^ 

B,{<P) = {ae Hi, \\a - 011^, ^ p, llall^, < A,] (2.25) 

where Ai is a sufficiently large constant. We also define Bp{M.i) = U^e;Kji3p(0). We 
choose this kind of neighborhood because in the following we need to control the norm in 
the uniform norm and in the a- Holder norm (for a < 

We now state our main result describing the dependence in e of the mean of the hitting 
time of a union of balls around the points of Al; starting from (pi^. 



Theorem 2.6. Under the assu'mptions \2. 1\ 2.4. for any minimum (pi^, and a set of minima 



Ail with Iq > I, there exists po such that for any po > p > 



where the error term satisfies "^{e) = 0{y/e \\n[e 



|3/2> 
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For the simple case where we have only three stationary points, two minima and one 
saddle (e.g. the case of the double well potential (11. 2p with Neumann boundary conditions, 
for 7 > I/tt^), we have the following corollary. 

Corollary 2.7. Let 0+ and 0~ be the two minima with S{(f)~) ^ S{(f)~^) and a the unique 
saddle point. There exists po such that for any po > P > 



- + ne)) (2.27) 

where the error term is '^{e) = 0{y/e \\n{e)f^'^) . 

2.4. Sketch of proof of Theorem 12.61 We first introduce the discretization we consider. 
The finite dimensional approximation of the SPDE is constructed as in the work of Funaki 
[27] and the work of Gyongy[28j. The approximation is defined via a spatial finite difference 
approximation of Equation (II. ip . 

We denote by the discretized potential, for y G R^"*"^ 

N 

SN{y) = hNYl wlrivi+i - y^f + y^Vi) (2-28) 

where /iat > is the step of discretization. We set = uo{xi) where uq G C;,c([0, 1]) is 
the initial condition and the Xi are the discretization points on [0,1]. Let us denote by 
Xi-i/2 the middle point of We construct a A^-dimensional Brownian motion B 

from the white noise W. Doing so we will be able to prove the convergence of to u in 
and almost surely. Thus we define, for 1 ^ z ^ 

Bl = -^W ( [xi_i/2, X [0, t]) . (2.29) 

The properties of the white noise imply that (-B*) are independent Brownian motions. 
The A^-dimensional process {Xt)t is the solution of 

dXi = -^VSNiXtYdt + J ^dBi for t = 1..N. (2.30) 

X^ and X^+^ are defined by the boundary conditions 

• for Dirichlet boundary conditions: 

X° = Xf +^ = 0, Vt ^ (2.31) 

• for Neumann boundary conditions: 

X° = Xl and Xf +1 = Xf , Vt ^ 0. (2.32) 

The discretized system is the linear interpolation between the points (xi,X*). To 
simplify, it is easier to adapt the parameters to the boundary conditions. 

• For Dirichlet boundary conditions, we choose 

^^=iVTl' ""'^nTi' VO^z^X + 1. (2.33) 
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For Neumann boundary conditions, we choose 
1 i 1 



h 



N 



Xi 



VO ^ i ^ iV + 1. 



AT' AT 2N' 

We set T^{B) the hitting time of a set B for the discretized system 

(fi) = inf {t > 0, n^(iV"H) e B} . 



(2.34) 



(2.35) 



We decompose the proof of Theorem 12.61 in several steps: 



1^1) for a given e and a sequence of initial conditions (pj^, each being a minimum 
of , converging to (pi^ (see Proposition 15. 6p . we prove that the expectation of 
{Bp{Mi)) converges to the expectation of the hitting time for the SPDE: 



W— 5>oo '0 " 



(2.36) 



To this aim, we use the convergence of u to the solution u. This is done in Section 

El 

(2) For a fixed A^, we compute the asymptotics of the transition time uniformly on the 
dimension. We get a prefactor a^i^) such that 



E^^[rf(Sp(M))] 



^P{e,N)<^{e) = 0{V^\\n{e)f') 



(2.37) 



where the error term "^{e) does not depend on N. This step is the main estimate 
and is detailed below. 
(3) The limit N ^ oo of aAr(e) gives us the correct asymptotics for the transition time 
in the infinite dimensional case: 



lim aAr(e). 



(2.38) 



This is done in Section |5l 

The estimate fl2.37p is proved in two steps. 

(i) First we start from a probability measure (the equilibrium probability: u^) on the 
boundary of a chosen neighborhood of the minimum (p^, which allows us to do the 



computation of aisf{e): 
1 



-E,M[Te{Bp{Mo))]-l 



(2.39) 



aN{€) 

This is done in Section O 
(ii) Then we have to control the error made by starting on the boundary of the minimum 
and not precisely at the minimum: 



aN{e) 



r,(S,(A^o))] - = < ^2(e) (2.40) 



with \&2(£^) = 0{^/e \\n{e)f^'^) . This result comes from the loss of memory of the 
initial condition adapted from Martinelli in [31] . This is exposed in Section |H 
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3. Discretization 



In this section, we present the convergence of the discretization to the solution of 
the SPDE and prove the convergence of the hitting times. 

3.1. Finite Dimensional Model. We write the discretized system in a mild form. 
We define a function kat, with [xj the integer part of x, 



[[N + l)x + \\ 

N +1 
[Nx\ + 1 _ 1 
N 2iV' 



for Dirichlet boundary conditions, (3.1) 
for Neumann boundary conditions. (3.2) 



We define the semi-group associated with the discretized Laplacian. The discretized 
Laplacian is a dimensional matrix, denoted by for Dirichlet boundary conditions 
and by for Neumann boundary conditions: 

/-2 1 



- — 







1 -2 

••. 









1 

-2/ 



A 



N 



"-N 



/-I 1 

1 -2 

••. 

Vo ... 













1 



(3.3) 



We consider the matrix p'^it) = hj^e'^"'^'^ . Therefore {t)ij is the solution of 

d 

P^(0),: 



(7A>^(t)),, 
1 



(3.4) 



The semi-group is the linear interpolation of p^{t) on [0, 1] x [0, 1] along the discretiza- 
tion points. 

Let us now prove the convergence of the solution of (13.81) to the solution of Equation 

Theorem 3.1. For all initial condition Uq G C^^([0, 1]), T > 0, and p ^ 1, we get the 

convergence 



.N 



U 



> u on [0, 1] X [0, T] 



(3.5) 



in the following senses: 

• m LP(fi,C([0, 1] X [0,r])), I.e. E 



I N 

\u — U 



VP 

loo,T 



^ 



almost surely in C([0, 1] x [0,T]), i.e. for every r] g]0, there exists E a random 
variable almost surely finite such that 
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Remark 3. Let us denote 



\u\ 



sup 

te[o,T] 



\u{x, t)f dx 



sup ||n(t)||^, . 

te[o,T] 



(3.7) 



gT ^ ii"'iioo T- As a consequence we get convergence in Theorem 13.11 in the L'^ 



We have IImIL -r. ^ \\u 
norm instead of the uniform norm. 

The convergence of the finite discretization is proved in [2S] if V is globally Lipschitz. We 
proved that the result holds in the case that V satisfies (12.51) via a localization argument. 
The idea, notably used by Funaki in [27], is to rewrite the finite dimensional system 
in a "mild form" and prove the convergence of this finite dimensional mild form to the 
infinite dimensional mild form (12. 3p . 

Lemma 3.2. For every uq G Cfec([0, 1]) and N > 0, the function defined on [0, 1] x M"*" 
satisfies the equation 



.N , 



U 



x,t) 



9t {x,K,N{y))uo{KN{y))dy 



Jo 



gt^^{x,KN{y))V'{u {KN{y),s))dyds 



+ V2ef [\ 

Jo JO 



g^^,{x,KM)W{dy,ds). {3.t 



For all p ^ 1 and T > 0, we have 



sup E 



TV 



sup 

[0,T]x[0,l] 



k^(x,t)r ^c(T,p). 



(3.9) 



Proof. This lemma is just a reformulation of the system of stochastic differential equations. 
We use the variation of the constant to integrate the linear part and then interpolate 
linearly the system to obtain a mild formulation of the function (see [271 EH])- To 
obtain the uniform moment bound, we proceed classically using a truncation procedure. 
We define and ur solutions of equations (13. 8p and (12. 3p in which we have replaced the 
function V by defined, for i? > by 

bniu) = V'iu)l[^R,n] + V'iR)l]n,+o.[ + V^'(-i?)l]-oo,-i?[. (3.10) 

bji is continuous, bounded and globally Lipschitz. Firstly, using the uniform estimates of 
the semi-group and the boundedness of bf>, we prove that for all T, all p > 1, there exists 
C {p, T, R) independent of such that 



sup E [|ti;^(x, t)!"] ^ C{p,T,R) < +00. 

[0,l]x[0,T] 

Secondly, there exists C{p,T,R) independent of A^, such that 



sup E 

N 



sup 

l-[0,l]x[0,T] 



U 



'x,t)f] ^ C{p,T,R) < +00. 



(3.11) 



(3.12) 



We use regularity of the solution (Kolmogorov's theorem) to prove (I3.12p . Thirdly, we use 
a comparison theorem to obtain uniform bounds on from bounds on u^^ where Rq is 
fixed and sufficiently large. □ 
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We use the convergence of to ur proved in [28] . 

Proposition 3.3. (^] For all R > 0, T > and < r] < ^ and uq in CfjO, 1], there 
exists a random variable almost surely finite such that 



I N I 



Proof of Theorem \3.1\ Let i? > 0, we define the stopping times 

tr = inf{t, \\uR{t)\\^ >R} = inf{t, 3x G [0, 1], \uR{x,t)\ > R} 
= inf{t, |K(t)|L > ^} = inf{t, 3x G [0, 1], > R}. 

Let us choose < 5 < 1. For i? > 1, we define 



(3.13) 

(3.14) 
(3.15) 

(3.16) 



First we show that P[f2ij] 
Proposition 13.31 for sufficiently large, 



{tr-5 > T and liminf Tn. > T}. 

TV— >oo 

> 1. Let M > 0. For a; G {^i? < M} n {r^j^^ ^ T}, by 



\u 



N\ 



(w) < \\uR\\^rj.{uj) + 6 < R 



(3.17) 



''i?lloo,r' 

which means that liminf Ar_s.oo t^(i^) ^ ^- Then by taking the complement relatively to 
< M} we get 

P[liminf <T;U< M] ^ P[r^_5 < T;^r < M] ^ r[rR^s < T]. (3.18) 

TV— >oo 

By definition of the time Tr^s, we have by the Markov inequality for p > 1 and from 
Equation (E 



P[liminf <T;U<M]^ F[tr.s ^ T] ^ P 

TV— >c« 



E 



\U\ 



oo,T 



{R - 6)P 



(3.19) 



Finally we get 



F[n%] = F[tr^s ^ T or liminf rjy ^ T] 

TV— >oo 

^ FItr^s ^T]+ P[lim inf <T;U<M]+ F[^r ^ M] 

N^oo 



2E[||n||^J 



(3.20) 



{R - by 

Since ^r is finite almost surely, we take first the limit M — )■ +oo then R — +oo. 

Let us define Q.r = Qr fl {^r < oo}. Since tr and are increasing in i? G N, the sets 
Qr are also increasing in R. Then we have 

nu^^.^R] = FlUR^n^R] = lim F[nR] = 1. (3.21) 

_R— >oo 

Let u G Qr. By definition of , there exists No{u) such that for all ^ A'^o(i^); 
(u) > T and tr-s{u) > T. By using the proposition 13. 3[ for all ^ No{u), 



oo,T 



\U 



'^-ur\\^^{uj)^U^)N-\ 



(3.22) 
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We define ^^w) by 



e^M= sup iv^lK 

Cji{uj) is finite on (Ir and is such that 
variable H by 



ML,Ti^)+U^)- (3-23) 
u\\ rr. ^ Let us define the random 



^^(w) on Qr \ Qr^i for i? ^ 2 
Cii^) on Qi. 



Then on Un^idji, set of probability 1, S is almost surely finite and \\u^ — u| 
which finishes the proof of the almost sure convergence. 



oo,T 



(3.24) 



To conclude, we show that E 
moments in (Lemma 13. 2p . we define 



N 

u — u\ 



IP 

loo,T 



converges to 0. Since IIm 



N\ 



loo,r 



R,No 



n 



N^No 



{tr_s > T and > T} 
We have = UjVofi/j.ATo- For all N Nq, we get by definition 



p 



1( 



■R.Nq 



I I 
\Ur - Ur\ 



oc,T + ^^kNo 







\K- 


- u\\ 



p 

oo,T ' 



Thus using Cauchy-Schwarz inequality and the bound (13. 9p . we get 



E 



IP 

,N 



^ E 



.N 



Ur 



IP 

loo,T 



Using the convergence of Ur to u (Proposition 13. 3p . we obtain 



lim sup E 

7V-s>oo 



p 

oo,T 



has uniform 
(3.25) 

(3.26) 
(3.27) 
(3.28) 



Let us fix ?7 > 0. Since P[f2ij] tends to 1 and VLr is increasing, we choose R such that 
P[f2^] ^ Tj. Similarly, ^r,No is increasing in iVo, thus P[f^ji;] = limAr^-^oo Pi^R atq] ^ 
Let us choose Nq such that P[f^K jv^] ^ 2?]. Inserting this bound in (I3.28p . we obtain the 
result. □ 

3.2. Convergence of the Transition Times. We conclude this section by proving the 
convergence of the transition times. 

Let us denote by uq the initial condition of the solution of Equation (11. ip and a 
continuous function. We define the hitting times: for p > 



r,(p) = inf{t>0,||n(t)- 011^ <p} 
rf(p) = inf{t>0,||n^(t)- 0^11^ <p} 



(3.29) 
(3.30) 



where 0^ is the linear approximation of 0. 



Proposition 3.4. Suppose that 
that for every pq > p > 0, 



lN 



6||^ converges to and that there exists po such 
E,Jr,(p)] < oo. (3.31) 
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Then for almost every p > 0, 



rf (p) -— ^ r,(p) a.s. and E,^[rf (p)] -— ^ E„„[r,(p)]. (3.32) 



N^oo N^oo 



Proof. For the sake of simplicity we omit e in the proof. First we prove that for all 5 > 0, 
T > 0, we have 

r(p + 5) A T ^ liminf r^(p) AT ^ limsup r^(p) AT ^r^p- 5) AT a.s. (3.33) 



From Theorem 13. 1[ \\u — u\\^^ converges to almost surely. Therefore with proba- 
bility 1, there exists Nq{u) such that for all ^ -^^0(0;) 

sup ||«^(t)-«(t)||^(a;)<^ and ||0^ - 0|| < ^. (3.34) 

Then for t ^ r(p + 6) AT and ^ A'o(a;), using the triangle inequality we get 

p + \\u{t) - 0IL ^ \\u{t) - t.^(t)||^ + 11..^ (t) - 0^11^ + \\u^ - 4^ 

<5+||«^(t)-0^||^ (3.35) 

which means that t ^ {p) A T . Thus, we obtain r(p + 5) A T ^ liminfAr^oo[T^(p) A T] 
almost surely. By the same arg uments for t ^T^{p) AT and N ^ Nq (cj), we get 

p^||n^(t)-0^|L^5+||n(t)-0|L. (3.36) 

Therefore limsup^^^Q [r^(p) AT] ^ r{p — 6) AT which proves the inequality fl3.29p . 

From the definitions of r(p) and t'^{p), the functions p H- r(p) and p H- T^{p) are left 
continuous and have right limits. Then using the fact that r(p) is finite almost surely, we 
get 

r(p+) ^ liminf r^(p) ^ lim sup r^(p) ^ r(p) < +00 a.s. (3.37) 

where r(p+) = lim^_5.o+ r(p + S). 

At a point of continuity of p H- t(p), we obtain r(p) = lim7v->oo t^(p)- Let us fix pi > 0. 
There exists ^ G Q a. null set such that for u ^ ^/K, p i— t- r(p)(u;) is bounded, decreasing, 
left continuous on [pi, +00 [. We define the set of discontinuities, 

^ = {(uj,p) eJ^'x [pi,+oo[,r(p+)(a;) ^ r(p)(a;)} C 1] x R. (3.38) 

Then we consider the projection from f2 x R on M along {u} x R. For u G we 
define 

I)(a;) = n^(^) = {p G [pi,+oo[,r(p+)(a;) ^ r(p)(a;)} C M. (3.39) 

P(a;) is at most countable since p t— )■ r(p)(ci;) is a bounded decreasing function. 

We define A/'(p) = 11^ (^) with 11^ the projection from x R on along Q x {p}. A/'(p) 
is the set of Q for which r(p) is not continuous at p. Therefore, we have 

^ = U^en {u;} X V{uj) = Up>p,Ar(p) X {p} . (3.40) 

Then, using Fubini-Tonelli Theorem 

/ P[Ar(p)]dp= / / l^(a;,p)dpdP(a;) = / / li,(^)(p)dpdP(a;) = 0. (3.41) 

J pi J Q, J pi J CI J pi 
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We get a null set S{pi) on [pi,+oo[ such that P[A/'(p)] = for all p G S{pi) i.e. the 
convergence is almost sure. To conclude, we consider a sequence (pn)n>o converging to 0, 
then S = Un^oS{pn) is a null set of M on which the convergence is almost sure. 

By using dominated convergence, we obtain the convergence of the expectations. □ 

4. Initial condition 

4.1. Large Deviation Control. For < a < 1, we set C"([0, 1]) the set of a-Holder 
continuous functions on [0, 1] equipped with the norm 

ll/llc^ = ll/lloo + sup '^f (4.1) 

x,y p y\ 

We also define -D"([0, 1]) the separable subset of this Holder space which is the closure of 
C°° in C". 

Let < a < I and p > 0, we consider the neighborhood -B" (</>) of </> G ^"^([0, 1]) 

B'^i<P) = {^e Dli[0, 1]), 110 - V^ll^. < p} . (4.2) 

We also have B^{Mi) = U^^m^B"^ {<!>)■ 

With this large deviation principle, Chenal and Millet [T7] derive exponential asymptotic 
estimates for the exit time of domains with a unique stable stationary point. Using their 
evaluations and the procedure developed by Freidlin-Wentzell [23] in the finite dimensional 
case, we have the following result. 

Lemma 4.1 ([IZ]). For < a < ^, there exists po such that for all p < po, we have for 
all (/) G B'^{(j)i^) and t] > 

limP^ [exp (^E-\S + 7/)) > TeiB'^iMi)) > exp (^e~\S - v))] = 1, (4-3) 

where S = S{(f)i^,Mi). Let n = Te{B'^{Mi)). Then 

where £ is an exponential variable of parameter 1. Moreover for all </> G Bp{(j)i^) 

lime log E,^ [r^] = S and lime log [r|] = 2^. (4.5) 

These estimates are the infinite dimensional version of the Freidlin-Wentzell theory. 

4.2. Exponential Contractivity. For a given tp^ = [tpi,--- jI^n) £ , we consider 
equivalently the point in and the function in C([0,1]) obtained by the linear inter- 
polation between the points {xi,ipi). Reciprocally, for ip G Cbc([0, 1]), we let ip^ be the 
linear interpolation of ip along the discretization. is the linear interpolation between 
the points {xi,ip{xi)). 

We set 

B^{^) = {^e CU[0, 1]), 11^ - 0IIOO < P} • (4.6) 
We adapt trajectorial results of contractivity for the localized process from Martinelli and 
Scoppola [51]. We denote u{(p),uji{(f)) the solutions of Equation fll.ip with respectively V 
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and bji, starting from 0. Accordingly, we denote u {(f) ),u^{(f) ) the solutions of Equation 
(13. 8 P with V and b^, starting from 



Lemma 4.2. Let (p be a minimum of S and R ^ Rq. There exists m,CR > and 
Eq, po > 0, such that for all p < po and every ip G B'^{(f)) we have, for all Eq > e > 



P 



sup 

N^No 



^ e 



-mt 



IIV^-0lloo>vt>o 



^ 1 



(4.7) 



This result can be proved via an adaptation of the arguments of [35] and [51]. Lemma 
14.21 describes that the solutions of Equation (11.11) and (13.81) depend slightly on the initial 
condition. Moreover, the solutions starting from two functions are exponentially close 
uniformly in the dimension. Martinelli and Scoppola called that the loss of memory of the 
initial condition because the specific initial condition is not relevant for the evolution of 
the process. 

4.3. Uniformity in the initial condition. Let us recall that is a minimum and M.i 
is a set of lower minima. We denote 



N I 



rf (5^" (0,, ) ) = inf { t , (t ) e 5," (0 J } 

(A^O = rf (i?^"(A^O) = inf {t, n^(t) G B';{Mi)] . (4.8) 

Similarly, we denote by rj^'^ the hitting time associated with the localized process u^. 

Proposition 4.3. For all pq > p > 0, there exists rj > such that for a sequence (pj^ of 
minima of , converging to in L? , 

s-v 



sup sup 



(4.9) 



For any sequence (f)f G of minima of converging to (pi G in , we also have 



sup sup 



(4.10) 



The proof comes from a comparison between the deterministic process (i.e. e = 0) and 
the stochastic process starting from the moment of the hitting time . 



Proof. Since the minima are not degenerate, we can assume p small enough to get 



5S 



^ -bU- 



(4.11) 



L2 



for some 6 > 0, all 1 < i < /, and all G B2p{(pi). 

First, let us prove similar estimates on the expectations of transition times for the 
localized process u^. We denote by a'^{(p^) the hitting time t^'^{A4i) for the process 
starting from 0^. We set 



Hr = { sup sup 



(4.12) 
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From Proposition W?2[ we get > 1 - e"'^«/^. 

Let us fix 6i > 0. We define T{e) = and we take e < Eq such that e~^'^'^^^ < p. On 
the set {a^(0jo) > T{e)], setting i) = we get 

(4.13) 



mT{e) 



with probabihty at least 1 — e Let us suppose that cr^ {(()) — a^{(f)if^) ^ and that 

The deterministic process m^' is the solution of (13.81) for the drift bf> and 5 = 0. (pf is 
a minimum of , so (pf is an equilibrium point of Then using Equation (14. lip , we 
get for t ^ 



L2 



^ e 



(4.14) 



L2 



by the triangle inequality. For t > to = ^ ln(16), we obtain 

From the large deviation principle, we can compare the deterministic solution with the 
perturbed one. We obtain C > such that 



P 



Ur (^A )-UR{ij ) 



P 



< ^ 
oo,2to O 



> 1 - e-^/^ 



(4.15) 

^ 5p 

r 9 ^ a 



Therefore, with probability at least 1 — e ""/^ — e we get ||u^(^/^)(2to 

which implies 

(a^(0)-a^(0j)+^2to. (4.16) 
We proceed similarly if o"^(0) — a'^^cpi^) ^ 0. In this case, we stop the process at a^{(f>). 
Finally we get \a^{(j)) — cr^(0«o)| ^ 2to with probability at least 1 — e~^'/^, for some C > 0. 
We obtain 

E[|a^(0)-a^(</)J|] ^E[|a^(0) 



cr 



Ni 



lo) \ lf^fll{a'V(0,^)>T(e)} 



+ E [|cr^^(0) - a^^(0j| (1^.^ + 1|..(^,^)>^(,)}^) 
^ 2to(l - e-^'/^)P [fin, n {a^(0j > r(£)}] 



(4.17) 



E 



|a^(0)-a^(0,,)| 

By using Proposition 14. 2[ we have P[f^^] < e"*"^/^. From Proposition 14. H we deduce that 
for e ^ Eo 

P[a^(0,J ^ Tie)] < 1 - e-'^"^ < e-^. (4.18) 
Moreover, we have for all ^2 > 



E 



a^(0)-a^(0j 



Ni 



(4.19) 



So we finally get 



E[a^^(0)-a^(0j]^2to(l-e- 



e-^) + e^(e-^/2^ + e'^) ^ e^. (4.20) 



By choosing 6i, 62 and r] small enough, we prove the proposition for the localized process. 
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Let US now choose R such that S{B'^{0), B'^{(f)i^^)) > 5* + 1, then from Proposition 14. ![ 
we have 

sup P<^[r,(5^(0)) ^ exp((5 + 1 - 6s) /e) = T^ie)] ^ e"^/^ (4.21) 

sup P<^[rf (A^O ^ n{e)] ^ e^^/^ (4.22) 

We consider the process u starting from </> and (pi^. Before T2{e), with high probabihty, 
the processes are in 5^(0) and coincide with up to this time. Moreover T2{e) is much 
larger than the transition time, so the transition aheady occurs when the processes reach 
B'^{OY. Therefore, with very high probabihty, the transition time for the locahzed process 
is exactly the correct transition time. 

For Equation fl4.10p . we follow a similar method, by using Proposition 14.21 for the local- 
ized process and then comparing the deterministic and stochastic processes in the neigh- 
borhood of a minimum. □ 



5. Approximation of the potential 

In this section, we prove (or refer to) results about the convergence of the potential and 
its related quantities. 



5.1. Convergence of the potential. Let us recall from Section 14.21 that for a point 
6 M^, we denote also by the linear interpolation between the points {xi,uf). For 
a function u G Cf,c([0, 1]), we denote by the linear interpolation between the points 
{xi,u{xi)). We say that the sequence G converges to m G if the sequence of 
linear interpolations associated to (also denoted u^) converges to u in the norm. 

Let us recall that HS^{u'^) is the Hessian matrix of at and can be interpreted 
as a bilinear form. We prove the following proposition. 

Proposition 5.1. For any sequence G converging to u & , we have 

• S^(m^) > S{u) < oo 

• for any sequence converging to h: VS^{u^) ■ y DuS{h) 

• for any sequences /i^, converging to h, k: 

ff5^(M^)(/i^, k"") y DlSih, k). 

N-^oo 

If u is twice differentiable DuS{h) = ^{u)h and if k is twice differentiable DlS{h, k) = 
Jq hHuSk. 

Proof. Let G be a sequence converging to n G , then converges uniformly on 
[0, 1] to u, so by dominated convergence. 
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The convergence in directly ensures us that 

TV X I 

^Y.^\<i-nff= I {u^)'{x)\x--^ f \u'{x)\'dx. (5.2) 
Let /i^ G be some sequence converging to /i G then we have 

i=l * i=l 

> / -fuh' + V'{u)h (5.3) 





by convergence of the derivatives and dominated convergence. Lastly, the convergence 
of the Hessian is completely similar. □ 

5.2. Convergence of the eigenvalues. Let us consider a sequence of points G 
converging to u in H^. We need to estimate the convergence of the eigenvalues 
{NXk,N)i^k!iN of ■ HS'^{u^) to the eigenvalues (Afc)i<fe of HuS. 

The convergence of a single eigenvalue NX^^j^ for k fixed, is obvious from Proposition 
15.11 The control of the convergence for all the eigenvalues is complex because of the higher 
eigenvalues (e.g. Xn,n)- This problem is closely related to the discrepancy between the 
eigenvalues of ;^A^ and 7 A, the discrete Laplacian (defined by f l3.3p ) and the Lapla- 
cian. We denote A^^, A° their respective eigenvalues in the increasing order. For Dirichlet 
boundary conditions, we have 



ek,N = NX%, - A° = 7 



(5.4) 



2N J 

Then eAr,TV = 7A^^(4 — tt^) does not converge to 0. The following proposition adapted from 
[30] gives us a control of the approximation of the eigenvalues and eigenvectors. 

Proposition 5.2. Let us consider a sequence mat G converging to u & and such 
that \\un " "wIIoo ~ ^ iw) ■ have: 

(i) there exist a G [0, 1[ and a constant Ci such that for all N and k < aN 

|A^AAr,fc - Afc - efe,Ar| ^ — , (5.5) 

(ii) there exists a constant C2 such that \eN^k\ ^ C2k^N~'^, 

(iii) for a fixed k ^ A^, the normalized (inH^) eigenvector (pk^f^f of HS^ {u'^) associated 
to Xk^N converges in to the eigenvector (p^ ofH^S associated to X^ and we have, 
for all k 



||0fc,7v||2,^ VN' 

Proof. The proposition is an adaptation of the results of [30] in our case since NHS'^ {u^) 
is the finite difference approximation of the Sturm-Liouville operator TiuS. The original 
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statement in [30] concerns an approximating sequence u which is precisely the sequence 



u 



N 



of hnear interpolations of u. If we take a sequence m , then for all y E 



N 



N 



U\ 



I l|2 



i=l 



(5.7) 



Since Hmtv — m||oo = ^iw)' deduce that the difference between the eigenvalues of 
N HS^ {u'^) and NHS^ {u^) is bounded by 0(-^) which gives us the result. A similar 
control holds for the convergence of the eigenvectors. The last result (15. 6 p comes from the 
fact that for the eigenvectors of T-LuS ([E] pp. 334-335), we have a constant C such that 
ll'/'fclloo ^ ^ \\'Pk\\L2- Then, since (f)k,N converges in H^, it converges in L°° and L^, then the 



result comes from the fact that ||0fc,Ar||2 7v ^ C^/N 



Remark 4. The normalized eigenvector e^r 



oc,N 



\\4>n\\: 



'>k,N\\i2- 

satisfies 



c 



N 



^n\\12 " N- 



□ 

(5.8) 
for the eu- 



The following proposition from [30] states uniform estimates in the function (p of the 
eigenvalues of the Hessian operators Ti^S and HS^{(f)^). 

Proposition 5.3. Let (pf, 0^ be sequences converging in to 0i, 02; then for all N, k 



Un\\2,N ^ 

Thus, this proves that the coordinates of the normalized eigenvectors in 
clidean norm are uniformly bounded by O (^"^ j ■ 



^ c 



I aJ. 



^ c 



(5.9) 



1,2. 



and XI = TT^P + V"{(j)i{x))dx + O (^) for i 

Remark 5. This proposition shows the convergence of the infinite product of the ratio of 
eigenvalues denoted by -D(0, ip) 



N 

n 



Afe(0) 



N 

n 

fc=i 



smce 



Afc(0)-Afc(V>) ' 
Afc(0)-Afc(V) 



Afc(V) 



C 

^ — • 



(5.10) 



(5.11) 



5.3. Product of eigenvalues. We show the convergence of the product ratio of the eigen- 
values of HS^{4>^) and HS^{^^) to D{(P,iP). 

Proposition 5.4. For any 0^, tp'^ converging in to 0, such that 'HS{ip) and 'HS{(p) 
do not have a zero eigenvalue, and that 

C 



V - i) 



loo ^ jY-2' 



(5.12) 
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we have the convergence 

det(i7^^(0^)) 



+00 



det{HS^{ilj^)) N^oo 



Afc(0) 



(5.13) 



Proof. The proof of the convergence comes from the fact that for small k the approximated 
eigenvalues are close to the continuous ones {Xk,N ~ Afc) whereas this is not the case for k 
close to (Proposition 15. 2p . The eigenvalues Xk,N{(p), ^k,N{4') are close at the first order 
in k uniformly on 0, ip (Proposition 15. 3p . Therefore we decompose the product in two 
parts for small k (i.e. k < aN from Proposition 15. 2p and large k. 

Let us denote Hk,N{<P) = NX]\f^k{(p^) — ^k{(p) — ^k,N- From Proposition 15. 2[ there exists 
< a < 1 such that, for k ^ aN , |/ifc,7v(0)| ^ 772- The same holds for the sequence ip^ . 
Then, we get, 

NXk,N{(l>) Afc(V') i + OkM) , , 0k,N{4>) - hAi^) 



where 6k,Ni<P) 



Afc(0) NXk,NW ^ + dk,NW l + OkM^) 

Xk{<p)~^{^k,N + fJ'k,N{(p))- Let us remark that for k ^ aN 



\Ok,Nm ^ 



C f k^ 



- 



A;2 \ iV2 



1 

iV2 



^Cla 



1 

iV2 



thus if we take a small enough and N large enough, we have |6'fc,Ar(V^)| < \. 
obtain 



aN 

mn 



NXk,N{<P) Afc(^) 



k=l 



Xk{<P) NXk,NW 



aN 



^2j2\0k,N{<i>)-0k,Nm 



k=l 



2Ca 
1^ 



(5.14) 

(5.15) 
Hence we 

(5.16) 



Ar2- 



since from Proposition 15. 3^ \0k,N{4>) ~ ^k,NW \ ^ 
For k > aN we proceed similarly. Let us write 

NXk,N{<P) Afc(^) ^ l + 0'k,N 
Afc(0) NXk,NW 1 + 



9[ 



k,N 



9i 



(5.17) 



where d'^^ ^ = Xk^Nii^) HAfc,Ar(0) — ^k,N{4')) and alike for 9[. From Proposition 15.31 we get 



for all k and N > Nq, that \9'k n\ ^ l^'^l ^ Thus we obtain 



fc2- 



mn 



iVAfc,^(0) Afc(V') 



Afc(0) iVAfc,^(V^) 



AT 



k=aN 



1 + T^ 



which finishes the proof. 

In fact, we need a slightly different convergence. 
Corollary 5.5. Let be 0^, tp'^ converging to (f), ip such that 

• N 



c 

< — 

N 



(5.1^ 



□ 



Then we have 



det (17^^(0^)) 



C_ 



(5.19) 
(5.20) 
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Proof. From the previous proposition, we get that 

det(i/5^(0^)) 



N(J,N\\ N^+oo 



where (j)^ (resp. ip'^) is the hnear interpolation of (p (resp. ip). So we prove 



D 



N 



det(i/5^(0^)) \det{HS^{(P'')) 



N 



-> 1 



where ^^(0) = >^k,Ni(p ) {K,n{(I> ) - Kn{(I> ))■ From the fact that 



lN 



(5.21) 



(5.22) 



\l2 ^ § we 



obtain 



_ AN 



L2 



^ j^. Then for all y G M^, we have 



N 



1 ^ 



i=l 



N 



/ \ I I 1 2 ^ 

(2;i) Myil ^ 



l4,Af 



Therefore we get that |Afc,Af(0^) — ^k,N{4>^) 
Then, we obtain 



L2 

jpj2- The same holds for ip. 



\y\\ 



2,N 



C 



1 

^ - 



k2 iV3/2 ^ p^iV 2 



(5.23) 



for sufficiently large. 
Thus we get 

N 



fc=l '^^^^ fc=l 



fc=l 



^-t(^) 

Then let us fix 77 > 0, we have 

. - 1 r 



kWN 



7]^ 



^==1 ■ k=r)y'N 

Therefore we get lim supjv_>.oo [-Dat]] = which proves the proposition. 



(5.24) 



(5.25) 



□ 



5.4. Approximated stationary points. The last property we need to check is that 
for each stationary point of S, there exists a unique sequence of stationary points of 
converging to this stationary point. Moreover, to ensure the limit of the ratio of 
eigenvalues, this convergence has to be fast enough (see Corollary 15. 5p . To this aim, we 
have the following proposition. 

Proposition 5.6. There exist C, Nq, such that for all N > Nq, there is for each minimum 
(resp. saddle point) 4> of S a minimum (resp. saddle point) cp^ of such that 

C 



(5.26) 
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where (p^ is the linear interpolation of (f). 



Proof. Since by Assumption 12.4^ there is a finite number of saddles and stationnary points 
tlien we only need to prove tlie proposition for a given saddle or minimum. Let be a 
minimum, we prove that there is sequence 0^ of minima of such that 



c 

< — • 

N 



The result fl5.26p follows from fl5.27p since we already have that 

C 













L2 




oo 



(5.27) 



(5.28) 



In order to prove fl5.27p . we use a fixed point theorem. Let us consider the ball B^^^ 
of radius in the |Hl2Ar norm where C is a constant we will fix later. We want to find 



j° G B^i^ such that V {(p^ + 2°) = 0. In that case we will have 

^ < — Wzf < — 



= (pN + z and 

(5.29) 



By a Taylor expansion of the gradient we have 
where gi is the remainder which can take the form 



-1 /^ScAf ^ 1 

g,{z)= I (l-t)__(0^ + tz)z2dt 



dz} 



N 



{l-t)V"'{<p., + tz,)z1dt. 



(5.30) 



(5.31) 



Then we have for all 2;, y G -B, 



and 



G 



9^^) - 9^[y)\ ^ ^ \A - yt\ ^ \zi - yi\ ■ (5.32) 



2Cn 



N 



Ar3/2 



Let us also remark that since is a stationary point for the potential S", thus we have 
— 70"(xj) + V\(p{xi)) = 0. Therefore we get 



V5^(0^). - - (-70"(x.) + V"{<P{x.))) 



1 



- 7iV'(0(a:.+i) - 2(P{x,) + 0(a;,_i)) - 70" (x,) ^ 



iV2 



(5.33) 



For sufficiently large HS^{(p^) is not degenerate then z'^ is solution of the fixed point 
equation 

z^ = HS^{^)-\-VS^{$^) - (7,(z°)) = F(z°). (5.34) 
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The (2, iV)-norm of HS 
absolute value). Then 



N(AN\-l 



is bounded by the inverse of the smallest eigenvalue (in 



\m\\iM^ 



HS 

1 



2,N 



^ C2N. For zeB, 



we get 



N 



2,N 



N 

Els 

i=l 



(j2 



(5.35) 



for C sufficiently small. Therefore F{B^^^) C B^^^. We also have for z,y E B^-,^^, 

F{y) - F{z) = HS^{r)-K-9^{y) + 9^{^))■ 
Then 

AT 



\F{y) - mwlN ^ c,N' E i-^^(^) + ^'(^)l' ^ 



i=l 



N 



2 II _ ||2 



Thus F is a contraction for sufficiently large. By the fixed point Theorem, there exists 



G B^i^ solution of z^ = F{z^) which proves Proposition 15.61 



□ 



6. Estimates 

6.1. Description. In this section, we compute uniformly in the dimension the expectation 
of the transition times. We proceed as in [2] and use the potential theory developed in [9] . 
Let us consider the A^-dimensional diffusion 



dYt = -VSN{Yt)dt + V2edBt (6.1) 
which comes from (I2.30p with the time change Y^^t = Xf. We denote by fi^ the invariant 



measure for the process Y 

yU^(dx) = e ^dx. 
Let us consider the norms for y G and p ^ 1 

N 



\y\\p,N 



x=l 



yi 



\\y\ 



00, N 



max \yi 



i=l-N 



(6.2) 



(6.3) 



Remark 6. As in the previous section, we associate to a point y G its linear interpolation 
on [0, 1] between the points {xi,yi) (xj is given by fl2.33p . fl2.34p ) that we denote by y. Let 
us consider the norm of y on [0, 1], we have for all p G [1, +00] 



(4A^)Vp ii^lip.^^ 



\y\ 



LP 



\y{x)f dx 



1 



jyl/p lll/|lp,Ar • 



This can be done using the Riesz-Thorin Theorem, remarking that 

and \\y\\oc.N=\\y\\i 



1 1 

y\\i,N ^ Ml^ ^ — \\y 



(6.4) 



(6.5) 
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In order to introduce the other norms, we need the following a priori estimates on the 
eigenvalues of the Hessian of S'^ . Let us recall the Hessian of at a point 0^ G is 

= -^(A^/i), + for h e (6.6) 

with the suitable boundary conditions. 

Lemma 6.1 ([SO])- For all (p^ e such that ||0^||^ < A, the eigenvalues {Xk,Ni<P^))k=i 
of HS^ {(j)^) arranged in increasing order satisfy the hound 

m{A)e - 1 ^ N\k,N{(p^) ^ M{A)k^ + 1 (6.7) 

where m{A) and M{A) do not depend on N and ip^ (only on A). 

Let us fix (p^ G M^. We consider the orthonormal eigenvectors {vi)i of HS^{(f)^). 
The decomposition oi h E in this orthonormal basis is given hj h = hiVi. For 
'p G [1, oo], we define the norms 

N 

Ml,T = Y.H \\h\\o.,T = ^^~h, . (6.8) 

i=l 

As in [2], these are the norms we use to control the approximations of the potential around 
our stationary points. Let us note that the norms depend on the point 0^. 

Remark 7. As in Section 4.1.1 in [2], the Hausdorff- Young Theorem can be adapted to the 



norms 



ip and II "lip all 2 ^ p ^ +oo and q such that q ^ + p ^ = 1, we obtain 



1|W|J^„«C(^|M|„T. (6,9) 



In fact, let T : — be the linear mapping T{y) = Yl!k=o Vk'^N By definition, 

\\Ty\\ p,T ~ II^/ILtv- '^^^ proof of (|6i9j) is an application of the Riesz-Thorin Theorem, 
between p = 2 and p = oo. On one hand, we have ||T?/||2 ^ = ||y||2 at since the eigenvectors 
form a orthonormal basis. On the other hand, we have HT'yllooAf ^ IIi/IIitv since the 
coordinates of the eigenvectors of the basis are bounded by (see Lemma 15. 2[ Equation 

O). 



Let us recall the infinite dimensional situation. The process u starts from a minimum 
bi^ of S and reaches the set of minima M.i. We denote by 5*0 = S{<pl^^,M.l) the height of 



the saddle points defined by (12.141) . 



we have a finite set hA^ = {x*} of 
we deduce that a sequence 
of such that each 



By Assumption 12. 4[ for all sufficiently large 
minima of S'^ . From Proposition 15.21 and Proposition 15.6 
of minima x'^^ converges to Similarly, there is a subset 
minimum of Aif converges to a minimum of A^;. 

We construct a graph for the finite dimensional the graph for the infinite dimen 

The vertices are the minima Ai^ 



sional case in Section 12.31 
points zl of for which So - S^{zl) 



The edges are the saddle 
< rj for some fixed 77 > 0. We connect the edge zl 
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between the two minima that the saddle point zl connects directly. To each saddle point 
zl, we associate a weight 



|Aiv(4)|e 



wi = ' (6.10) 

^\deiHS^zl)\ 

To each minima x*, we associate a value = a{x*) G M. We denote by aj+ and aj„ the 
two values associated to the minima connected by the saddle point z*. 

We associate to this graph a quadratic form Q^{a), for a a real vector indexed by the 
minima 

^'^(«) = 5Z«^r(«;+-ai-)'- (6.11) 

The equivalent conductance, C*{N,e), between the sets x^^ and Aif is defined by 

C*{N,e) = inf {g^(a),a(x;j = l,a{x*) = 0,x* G Mf"} . (6.12) 

We recall the fundamental formula (16.151) proved in [9]. The expression of the expectation 
of the hitting time t^{B^ {^Iq)) is based on two quantities: the equilibrium potential and 
the capacity with respect to the sets (x^^J and (J^f). The equilibrium potential, 
h*, is defined by h*{x) = P^[rf (S^(x;j) < r^^ {B^ (Mf))]. The Dirichlet form, (f^, 
associated with the diffusion process Y on is 



>g'\h)=e ||V/i(x)||2^/i^(dx). (6.13) 

The capacity is the evaluation of the Dirichlet form on h* . The capacity also satisfies a 
variational principle. We have 

= inf {(^^(/i), /i G ifi(M^), /i = 1 on B^{xl), h = on B^{M^)} . (6.14) 
The expectation of the hitting time is expressed by 

where is a probability measure on dB^ i^lo^- 

6.2. Capacity. We prove that the capacity defined in (16.141) can be estimated by the 
equivalent conductance C*{N,e) defined in (16.121) . 

Proposition 6.2. For all e < Eq and p, we have 

cap{B^{xl),B^iMf))=eV2^e " C*(iV, e)(l + ^i(e, iV)) (6.16) 

where lim supjv^_,_o^ \^lJi{e,N)\ < ^/e \\n{e)f^'^ for all N > Nq. 

The proof of this result is an adaptation to the case of a finite number of saddle points 
of Proposition 4.3 in [2]. The estimate of the capacity is made in two steps: an upper 
bound and a lower bound. 
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6.2.1. Upper bound. We have the following proposition. 
Proposition 6.3. For all e < Eq and p, we have 



ce^p{B^{xl),B^{Mf')) ^eV2^ ' C*{N,e){l + N)) (6.17) 
where limsupjy^o^ |?/^„(£:,iV)| < ^/e \\n{e)f^'^ . 

Proof. The proof of this upper bound follows the proof of Lemma 4.4 in [2]. To obtain an 
upper bound for the capacity, we just estimate the Dirichlet form on a test function h~^. 
h'^ is defined on some neighborhood Cf{z*) of each saddle point z* for some 5 > small 
enough. 

In the local orthonormal basis (given by coordinates i/^*-* G M^) of the saddle point 2;*, 
the neighborhood Cf{z*) is defined by 

Cf{z:) = L« G : ^ S-^, ^ / ^ iV - 1 1 + (6.18) 

3/2 

where (r;) is a sequence satifying < 00 and {\n,i)i are the eigenvalues in the 

increasing order of HS^{z*). Let us denote = UiC^{z*). 
Let us consider 

Sj^^s = {x, S'^ix) ^ S^{z*) + c5\\/i} . (6.19) 

The set {Sn,5 U CfY contains a finite number of connected components denoted Dj since 
each of them contains at least a minimum x* (which are in finite number by Assumption 
12. 4p . For each connected component Dj, we define to be the constant aj G [0, 1]. For 
a saddle z*, we denote Di+ and Di_ the connected components attained from z* when 
1/^*^ = ((5(Jo,0) and y^^^ = {—Sao,0) respectively. 

On Sn,5 \ Cg^ , we take h'^ of class and such that || V/i~''||2 at < ^- Then we define h'^ 

on each C^{z*) in the local coordinates, by h'^{y^'^'>) = fi{yQ^) where 
Therefore, we have to estimate S^{h^) = Ii{i) + /2 with 



h{t)=e / ||V/i+(x)||'^e-^dx, /2=W ||V/i+(x)||'^e-^dx. 



2,Af 

(6.21) 



Taking 6 = e |ln£|, the integrals /i(i) give us the right asymptotics and are estimated 
by an adaptation of Lemma 4.4 from [2]. The quadratic approximation of the potential on 
the sets Cf{zX) is a consequence of Remark [7] and of the choice of the sets C^{z*). The 
integral I2 is computed by following the same method as in Lemma 4.6 in [2]. 
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Therefore, we obtain that for all {cij)j, for ^ No{e) 



p V p V ' - ^ y/\det{HS^{z*))\ ^ iv i v ;i y 

Taking the minimum of the right-hand side over a, we get the result fl6.17p . □ 
6.2.2. Lower bound. We now prove the corresponding lower bound. 
Proposition 6.4. For all e < Eq and p, we have 

cap{l3^{x*),B^{M^)) ^eV2^ C*{N,e){l + ^i{e,N)) (6.22) 
where limsup^y^o^ \'ipi{e,N)\ < ^/e \\n{e)\^^^ . 

Proof. The proof is adapted from [2]. For a saddle point z*, we take a narrow corridor from 
one (local) minimum to another one and minimize the Dirichlet form on the union of these 
corridors. In [2] , this corridor was a rectangle because of the particular case considered. 
In this article, we have to be more precise about their construction. We use the same 
notations as in the proof of the upper bound. 
Let us fix 6o. We consider the subset of R^^^ 

^(^n = ly^'^ e : b«| ^ 5-^, 1 ^ / ^ iV- ll . (6.23) 

and we define C^{z*) = [-60,60] x C^''^{z*) + z*. We denote by x*_ and x*+ the two 
minima of the basins surrounding z*. 

Let (7o(s))sg[_s_,s+] be a regular path from Xj_ to Xj+ with 7o(s) = z* + (s,0) for 
s G [— (5o,5o]. We also suppose that there is r/ > for which 5^(70(3)) < Sq — 3// for 
|s| ^ 60 and that ||7o('S)||2 jy ~ 1- Let, for all s, A{s) be an isomorphism from R^~^ to 
70(5)-*- C R^ of class in s and such that for |s| < 60, A{s)y = {0,yi, ■ ■ ■ ,yN-i)- Then 
we construct a family of paths 'j{s,y±) by 

7(s,|/±) =7o(s) + A(s)yx. (6.24) 

Such a construction of a path 70 is always possible in the infinite dimensional setting 
(because of Assumption 12.41) . Then taking the finite dimensional projection, it gives us a 
path for the finite dimensional case. 

We define the corridor from Xj_ to Xj^, for 6 > small enough 

Cs{z*) = [x = 7(5,2/±),2/± G Cf'^«),V.} . (6.25) 

Let h be the equilibrium potential which realizes the minimum of the Dirichlet form and 
define a^^{y±) = h{xi± + A{±s±)y±), the values near the minimum. 

To estimate a lower bound, we are going to restrict the Dirichlet form on the union of 
the corridors Cs{z*): 

S'^{h)=e[ ||V/.||^,^/i^(dx)^^e / \\Vh\\l^f,'^{dx)=eJ2W. (6.26) 
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We define the function fi on Cs{z*), by fi{s,y±) = h{'y{s,y±)). The change of variable 
on Cs{z*) gives us the Jacobian gi{s,y±) = det{Jj){s,y±) and we obtain 

2 



5/. 



ds 



9i{s,y±)dsdy_i 



(6.27) 



We take y± as a parameter then the second term is bounded below by the minimum over 
functions /j of the integral 

2 

g-5'^(7(s,J/x))/e^.(5^^^)ds (6.28) 







/: 


ds 



with the conditions fi{—S-,y±) = h{xi- + A{—S-)y±) = a^~{y±) and fi{s+,y±) = h{xi+ + 
A{s-^.)y±) = a''~^{y±). This gives us a lower bound for the capacity. 

A simple computation shows that the function fi realizing this lower bound is 

Ms,y,) = ia^^y,) - a-iy,)) '^^- ,s^is,y,),eg^^^^y^y.^, + ^^-(^^)- (6-29) 
Inserting this function in the integral fl6.27p . we obtain 



h(i) > 



{a'^iy±)~a^-iy±)y 



dy^. (6.30) 



The end of the proof comes from an upper bound of the integral uniformly for y± G 



C^'^{z*). We write 



,s^(^'y^)/^g^(s,y^)-'ds = h{i) + m 



(6.31) 



where 

hit) 



r e'"^''''^^^'9^is,y^r'ds and hi{)=[ e^^^^'y^^/' g,is,y^)-'ds. (6.32) 

J -So J\s\>5o 

As in Lemma 4.8 in [2], we control the quadratic approximation near the saddle z* with 
the following lemma for which we omit the proof. 

Lemma 6.5. For all y = {s,y±) G C^{z*), if the sequence {ri)i satisfies < oo, we 

have for 6o ^ 6 



S''{l{s,y^) + z*) - ^^^(7(0,i/x) + z*) + - |Ao,;v| s 



2 



^^(.*+7(0,yx))-5^«)-^$^A 



N,kyl 



k=l 



< A,6\ 



(6.33) 
(6.34) 



Following the proof of Lemma 4.7 in [2], we can also prove the existence of a constant 
Aq such that for all and yj_ 



IS) ^ e 



lA 



1 + 



(6.35) 
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In addition, we need to prove an upper bound for the integral hii). 
Lemma 6.6. There exists a constant Aj such that for all N and y± 

Ijii) < AjVNe^ (6.36) 
where t] > is given by the definition of the path 70. 

Proof. We have to be careful with the change of variable. Let us write the Jacobian matrix 
Jl{s,y±) in the local base (7q(s), 7q(s)"'"), if we denote Pq the projection on Span{'jQ{s)), 
we get the Jacobian matrix (written by blocks) 

^0^(...x)^(^^^°<fW^-' (6.37) 

since Imyl(s) = 7q(s)-'". Then, as A{s) is an isometry, we obtain that 

g,{s,y^) = |det(J7(s,2/x))| = \1 + Po{A'{s)y^)\ = 1 + 0(5). (6.38) 



Thus, for 6 sufficiently small, 

/7(0 = / e^''(^'^^)/^(7,(s, y±)-'ds ^ (1 + CS)e^{s+ + s.) ^ 2{s+ + s_)e^ 

-'k|><5o 

since S^{s, y±) < S — 2ri for all |s| > 5q, and G C^'"*". Then by construction of the path 
we have that 

s+ + s„ ^ C - Xi+W^^^ ^ Ca//V - Xj+ll^a . (6.39) 

□ 

We insert f l6.35p and f l6.36p in Equation f l6.30p . Then we proceed as in the proof of 
Lemma 4.7 from [2] and we obtain 



l + A^ + A'c-'i 



(6.40) 



Using Equation fl4.10p from Proposition 14. 3[ we obtain for all y±, \a^{y±) — a^{0)\ < e~~. 
Then using the approximation fl6.34p and following the proof of Lemma 4.7 in [2], we 
obtain for 6 = \/ Ke |ln£:| and 5o = -^'^ |liie| with K' > K, 

-N~2 



/5(z)^e(a--a-+)V^^£ l^^'"' {l-A,^eMe)f\ (6.41) 
' v/|det(/75^«))r 5v I V ;i ; \ ) 

Equation f l6.22p follows by minimizing along the {a^)j- □ 
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6.3. Uniform estimate of the mass of the equihbrium potential. We prove esti- 
mates of the numerator of ( 16.151) . Let us denote G to be the closest minimum to 
(piQ in L'^{[0, 1]). We will prove an adaptation of Proposition 4.9 of [2]. 

Proposition 6.7. For all s < Eq and p, we have 

/i*(x)d/i^(x) = ^ ^ ^ =e r^(l + ^2(£,iV)) (6.42) 



'detHS^ixlj 
where \M^,N)\ < y/e\\n{e)f^ for all N > No. 

Proof. As the previous section, we define around the minimum x'^^ G a neighborhood 
Cj^(x;*^). In the local orthonormal basis of the minimum x^*^, the neighborhood Cj^(a;^^) is 
defined by 

C^K) = |y G : \m\ ^ ^^-/==^ ^ / ^ iV- l| (6.43) 

^3/2 

where (r;) is a sequence satifying J2i 7372 < 00 and {\n,i)i are the eigenvalues in the 
increasing order of HS'^lx] 
We need to estimate 



I W2 

lo)- 



/i*(x)d/i^(x). (6.44) 



Let us remark that for x G dC^{x*), then one of the coordinate is precisely Sr^/ ^jXk^N 
thus 

S^'ix) > S^ix*) + 6^1 - C6'' > S^ix*) + c6\ (6.45) 

We consider S' such that the set {0, S{(f>) E.]S{(j)ig), S']} contains no stationary point. Then 
using Proposition 15. 1[ for all 77 small enough, there exists A'o such that for N > Nq, 
[x, S^{x) G [S^{x*) + |c5^, S" — 77]} contains no stationary point. We define the set A = 
{S^{x) ^ S^{x*) + c6^}\ B^{x*). Note also that for S small enough, Cf (x*) C B^{x*). 
Hence we decompose (16.441) in three parts: 



h*{x)dfi^{x) = l8+ / h*{x)dfi^{x) + / h*{x)dn^{x) (6.46) 

To estimate the third integral we need a control on the equilibrium potential on the set A. 

Lemma 6.8. For all p < po and rj > there exists EqIp) such that for e < Eq and 6 > 0, 

let X E A, we have 

hUx) = P.[rf (i?;(x*)) < T^iB^iMm ^ e-(^'-^'^(^)-^'')/^ (6.47) 

Proof. By definition of the set A all the paths from x E A to x* attain a height of S' — rj 
at least. To prove this fact, let us take a path from x to x*, it must attain its maximum S 
at some time to. This maximum must satisfies 5* > ^^(x*) + C(5^, since if it is not the case 
then from Equation fl6.45p . the path must stay in C^{x*) which contradicts the fact that 
X is in A. Then the minimal path from x to x* must attain its maximum at a stationary 
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point of height greater than S'^{x*) + cS^ thus of height greater than S' — rj. This gives us 
an easy lower bound for the rate function on the set of transition from x & Aio x* . Then 
using the method from [2S] and the uniform large deviation principle, we prove that 

h*{x) = P,[rf (x*)) < r^{B''^{Mf))] ^ e"^^'"^'^-^"^^))/^ (6.48) 

uniformly in A^. □ 

We get from fl6.46l) 



ix)^S^ix* )+c<52 JS'^ix)!S,S^{x* )+c<52 



(6.49) 

where we have used the fact that h* is bounded by one for the second integral and the 
previous lemma for the third integral. The integral Is gives the main contribution and is 
estimated as in the proof of Proposition 4.9 of [2] using the quadratic approximation of 
the potential on C^{x*i^). The second integral on the right-hand side is estimated as in 
the proof of Lemma 4.6 in [2]. 

We bound the third integral by the volume of the set [^S^ {x) ^ S'^ {x\^) + c5^} which is 
bounded uniformly in N . In fact, from the bound on and the convergence of S^{x*i^) 
to S{(f)ig), we get for 6 sufficiently small 

{5^(x) ^ S'^ixl) + cS'}c { II V^x||^^^ + \\x\\l^ < iV(5(0j + c)} (6.50) 

which is a deformed ball. The computation shows that this quantity is uniformly bounded 
in A^. 

We obtain the result since the order of magnitude of the two last integrals (O (^e^(^'^'^y^^ ) 
of 06.491) is much smaller than 1^ = (^^~^ (^lo^/^^ □ 



6.4. Finite Dimensional Formula. The finite dimensional Formula is now obtained 

^0 



'^^ is the unique minimum of such that 



C 



where 0^ is the linear interpolation of 0;^. 

Proposition 6.9. Let he the transition time from {(pf^) to Bp{Aif), we have uni- 
formly in N 

V [^e] = , (1 + ^ie,N)) (6.52) 

C*{N,e)^detHS^) 

where C*{N,e) is the equivalent conductance and 

limsup|^(e,A^)| ^Cy^\\nsf\ (6.53) 
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Proof. Inserting the estimates for the capacity (Proposition |62]) and the numerator (Propo- 
sition 16. 7p in Equation fl6.15p we conclude that 



27re- 



C*{N,s)JdetHS^{<P^J 



(6.54) 



where hm sup^ |^i(e,iV)| < Cy^\\n{6)f^ and is a probabihty measure on dB^ (0z^)- 
Now we use Proposition 14.31 to replace the measure by the point For y G (0z^), 
we have by definition 



)-S^{y) <p. 



-y\\L2 <p 

Then from Proposition 15. 6[ we have Nq such that for N ^ Nq 

U,, ~ y\\% < 2p' |5(0,J-5^(y)| <2p. 



(6.55) 
(6.56) 



Thus since V is regular, we obtain 



,/ ||2 II ,||2 



Let z = y — (f)i^,, we have by integration by parts 





/l|2 1 


A' l|2 




\Wlo - 









2{<t>lz') + \\z'\ 



L2 



< Cp. 

-2{4>lz) + \\z'\\l, <Cp (6.57) 



since (pi^ is regular as a classical solution of a differential equation. Then we obtain by the 
Cauchy-Schwarz inequality 



Thus we get 



\z'\\% ^cp + 2 ii^y 1^, ^ (c + 2 ikii^op- 



y-<L<\\y~ ^doo ^c'\\y- = c \\z\\^, ^ c"^p. 



Using Proposition 14. 3[ we get that for all N ^ No 



^ e' 



S-2ri 



(6.58) 



(6.59) 



(6.60) 



which gives us (16.521) since the exponential asymptotics of fl6.54p is greater than e 



□ 



6.5. Proof of Theorem 12.61 From Proposition 16.91 applied to the finite diffusion approx- 
imation where the minima and saddle points are given by Proposition 15. 6[ we have 

V ['e] = ^"^f ^ (1 + ne,N)) (6.61) 

where the factor hj^ comes from the time change (Equation fl2.30p ). Using Proposition 15.21 
(convergence of the eigenvalues) and Corollary 15.51 (convergence of the ratio of eigenvalues) , 
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the quadratic forms converges to Q: 



g^(a) Jdetifg^«) -— ^ V \\^{<P\)\ J . , e-^(a,+ - a,.f 



S(<l>*i) 



Q{a)e —J-Dem^, S. (6.62) 



where are the relevant saddle points given by Proposition 15.61 Then the minimizer 
converges. For all e, we get 

1 I zr s(4.*,) 

lo . 



-C*(iV,e) Jdet//^^(0;^) -— ^ C*i<t>io,Mi)e-^jDetn^^^^S. (6.63) 

AItv iV— >oo 



Therefore, we obtain the result of Theorem 12.61 from Proposition 13. 4[ 
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